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Preface 

The  investigation  reported  herein  was  sponsored  hy  the  Defense  . 
Nuclear  Agency  (DNA)  under  Subtask  Y99QAXSB211,  "Penetration,"  Work 
Unit  l4,  "Technical  Penetration  Developments."  The  development  of  a 
rock  penetration  nomogram,  which  was  the  end  product  of  the  investiga- 
tion, was  suggested  by  MAJ  Todd  D.  Strong,  CE,  who  was  the  DNA  Project 
Officer  at  the  beginning  of  this  work  unit. 

The  study  was  conducted  by  personnel  of  the  Soils  and  Pavements 
Laboratory  (S&PL),  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  during  July-October  1976  under  the  general  supervision  of 
Messrs.  J.  P.  Sale  and  R.  G.  Ahlvin,  Chief  and  Assistant  Chief,  respec- 
tively. Mr.  R.  S.  Bernard  developed  the  penetration  formula  with  the 
technical  guidance  of  Drs.  P.  F.  Hadala  and  B.  Rohani  under  the  direct 
supervision  of  Dr.  J.  G.  Jackson,  Chief  of  the  Soil  Dynamics  Division, 
S&PL.  Mr.  Bernard  also  prepared  the  report. 

COL  J.  L.  Cojinon,  CE,  was  Commander  and  Director  of  WES  dm'ing  the 
investigation  and  at  the  time  of  preparation  of  this  report.  Mr,  F.  R. 
Brown  was  Teciinicol  Director. 
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EMPIRICAL  ANALYSIS  OF  PROJECTILE  PENETRATION  IN  ROCK 


Introduction 

The  design  and  deployment  of  earth  penetrating  weapons  requires 
simple,  reliable  techniques  for  predicting  their  performance  in  soil 
and  rock  targets.  Bernard  and  Creighton  (Reference  1)  have  developed 
a time-sharing  computer  code  (PENCO)  that  analyzes  projectile  motion 
within  layered  earth  targets.  Although  the  theory  used  in  PENCO  is 
fairly  simple,  the  calculations  are  too  cumbersome  to  be  performed  by 
hand.  It  is  therefore  desirable  to  have  an  even  simpler  prediction 
technique  that  does  not  require  a computer. 

Young's  equation  (Reference  2)  has  proved  to  be  the  simplest  and 
most  effective  empirical  formula  for  soil  penetration.  This  equation 
uses  an  empirical  penetrability  index,  which  varies  from  target  to  tar- 
get and  which  must  be  determined  from  previous  penetration  data.  Soil 
penetration  has  been  extensively  investigated  in  recent  years,  and  the 
inforaation  already  on  record  is  adequate  for  generating  reasonable 
pi'edictions  in  many  types  of  soil. 

Young's  equation  has  occasionally  been  used  to  correlate  rock 
penetration  data,  even  titough  the  equation  itself  was  derived  from  soil 
penetration  data.  Since  it  has  never  been  sUowt^  that  the  scaling  rela- 
tions for  soil  are  the  some  as  those  for  rock,  thero  is  clearly  a need 
for  on  equation  based  on  i^enotration  data  for  rock  in  situ.  The  de- 
velopment of  a penetration  fomula  specifically  for  rock  is  summarized 
in  the  following  paragraphs. 

Efttpirieal  Rock  Penetration  Equation 

The  first  stop  in  an  empirical  analysis  is  to  quantify  the  effects 
of  those  variables  that  c:ui  bo  specified  most  accurately.  Por  penetra- 
tion, these  variables  ai'e  the  pi'ojectile  mass  (M),  diameter  (D),  nose 
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shape  (CRH),^  and  impact  velocity  (V).  Once  the  relation  "between  pene- 
tration and  projectile  parameters  is  understood,  the  analysis  can  pro- 
ceed to  the  effects  of  variables  that  are  less  accurately  specified 
(i.e.  the  target  parameters). 

Experimental  data  for  concrete  (Reference  3)  indicate  that  the 
final  penetration  depth  (Z)  is  roughly  proportional  to  the  projectile 
mass  and  impact  velocity  and  inversely  proportional  to  the  square  of 
the  projectile  diameter.  Most  rocks  are  hard  ajid  brittle  (like  con- 
crete), so  it  is  assumed  that  the  same  trends  apply  for  intact  (un- 
jointed) rock.  Of  course,  the  degree  of  penetrability  will  change  from 
target  to  target. 

Penetrability  is  probably  related  to  a number  of  target  variables, 
such  as  density,  shear  strength,  tensile  strength,  bulk  ajad  shear 
moduli,  and  other  properties.  Unfortunately,  it  is  rare  to  find  pene- 
tration data  for  rock  sites  where  all  these  quantities  were  known  with 
consistent  accui'acy,  even  when  the  rock  was  in  an  intact  condition.  It 
is  therefore  assumed  that  the  penetrability  of  intact  rock  depends  only 
on  the  mass  density  (p)  and  the  unconfined  compressive  strength  (Y), 
which  are  the  two  properties  most  often  measured. 

Natural  rock  deposits  often  contain  joints  and  fractures,  so  that 
the  unconfined  compressive  strength  of  intact  core  specimens  may  fail 
to  characte^.'ize  the  rock  as  a whole.  Deoi’e  (Reference  1*)  has  proposed 
the  Rock  Quality  Designation  (RQO)  as  on  in'^ax  for  the  degi’ee  of  frac- 
turing of  the  rock  in  situ  at  a given  site.  'The  value  of  the  RQD  is 
determined  by  a modified  core  logging  procedure: 

All  solid  pieces  of  core  that  are  10  cm  long  or  longer  are 
added  up,  and  this  length  is  called  the  modified  core  re- 
covery. The  modified  core  recovery  Is  divided  by  the  total 
length  of  core  run,  and  the  quotient  multiplied  by  100  per- 
cent is  the  value  of  the  RQD. 


pmi  ogive  nose  radius 
**  projectile  diameter 


The  RQD  has  previously  been  shown  to  be  related  to  other  qualitative 
rock  descriptors,  such  as  fracture  frequency.  Now  it  seems  also  to  be 
ein  index  of  penetrability. 

A number  of  dimensionless  equations  can  be  constructed  that  relate 
the  projectile  and  target  variables.  One  of  these  has  the  form 

^ = f(V.^  , RQD  , CRH)  (1) 

3 

All  quantities  in  Equation  1 must  be  expressed  in  compatible  units.  The 
left-hand  side  represents  a dimensionless  penetration  depth,  and  the 
first  group  (V/p/Y)  on  the  right-hand  side  represents  a dimensionless 
impact  velocity.  Aside  from  the  introduction  of  the  RQD»  Equation  1 
contains  natural  dimensionless  groups  involving  the  mass,  velocity, 
diameter,  and  geometry  of  the  projectile,  as  well  as  the  intact  strength 
and  density  of  the  target. 

Physical  property  data  ore  incomplete  for  the  targets  in  most  of 
the  rock  penetration  tests  on  record.  However,  tai'get  descriptions  ade- 
quate for  the  present  analysis  do  exist  for  11  rock  penetration  tests 
conducted  by  Sandia  Laboratories  under  EBDA  or  DNA^  sponsorship  (Rofer- 
orces  ^-9).  Figure  1 shows  these  data  plotted  with  concrete  penetration 
data  taken  fi’om  Reference  3.  Table  1 enumerates  the  23  individual 
test  results,  and  Table  2 gives  the  projectile  and  target  parameters. 

The  combined  data  were  used  to  obtain  a penetration  equation  of  the  some 
general  form  as  Equation  1. 

The  penetration  data  for  concrete  (Coiifiold  and  Clator,  Refer- 
ence 3)  and  for  welded  tuff  (Sandia/DNA  tests.  Reference  $)  constitute 
Che  baseline  for  intact  rook  (RQD  “ lOO)  in  Figure  1.  The  other  data 

^ For  example,  p could  be  expressed  in  g/cm^,  Z ajid  D in  cm,  M in 

. grams,  A in  cm2,  y in  em/s,  ai\d  Y in  dynos/cm^. 

ERDA  is  an  ancrotiyw  for  Energy  Research  ojid  ^velopwent  ^oney, 

e and  DMA,  for  Defense  Nuclear  ^eney. 

■ The  rock  penoTration  data  are,  by  themselves,  a bit  too  scat- 
tered for  comfort.  The  addition  of  the  concrete  data  fills  in 
the  loose  pattertt  of  the  rock  data,  facilitating  a reasonable 
con'olation  of  all  the  test  results. 
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(for  RQD  < 100)  are  documented  in  References  5 and  6.  The  properties 
for  the  concrete,  the  welded  tuff,  and  the  82-RQD  sandstone  are  dis- 
cussed in  References  3,  7,  and  9»  respectively.  The  properties  of  the 
remaining  targets  are  given  in  Reference  9* 

Various  attempts  were  made  to  correlate  the  dimensionless  param- 
eters that  appear  in  Equation  1,  Figure  1 shows  the  most  successful 
normalization  of  the  data.  The  equation  for  the  dashed  line  in  Figure  1 
is 


M/A 


(2) 


7T  2 

where  A = is  the  projectile  cross-sectional  area. 

The  error  bars  in  the  figure  reflect  t^e  relative  uncertainty  in 
the  value  of  the  target  strength.  There  are  many  test  results  on  record 
that  could  not  be  used  in  the  formulation  of  Equation  2,  due  to  a lack 
of  descriptive  data  for  the  targets.  However,  the  points  in  Figure  1 
represent  tests  for  which  the  target  density,  strength,  and  RQD  were  all 
known  within  reasonable  bounds. 

The  available  data  did  not  show  a significant  increase  in  penetra- 
tion depth  as  CRH  increased  from  l.§  to  9*2^,  so  the  nose  shape  variable 
was  dropped  from  the  analysis.  However,  the  equation  is  not  considered 
valid  for  blunt-nosed  projectiles,  since  no  blunt  projectile  data  wore 
used  in  its  formulation. 

Expressed  in  metric  units.  Equation  2 becomes 


where 

Z a final  penetration  depth,  cm 
M « projectile  mass,  kg 
D projectile  diameter,  cm 

V “ Impact  velocity,  m/a 

p a target  bulk  density,  g/cm 

Y a unconfined  compressive  strength  of  the  intact  rock,  bars 
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RQD  = Rock  Quality  Designation,  pet 
Equation  3 is  expressed  as  a nomogram  in  Figure  2,  which  includes  a 
sample  calculation  using  the  nomogram. 

Limitations 


Although  concrete  data  were  used  in  thy  development  of  Equation  2, 
this  expression  is  not  intended  as  a concrete  penetration  formula.  The 
concrete  data  were  introduced  only  to  help  establish  the  baseline  for 
intact  rock.  There  are  several  equations  already  docximented  for  con- 
crete (e.g.  Reference  10),  which  produce  better  results  in  concrete  than 
Equation  2.  For  rock  targets,  the  penetration  data  are  more  scattered 
(and  the  target  parameters  more  uncertain)  than  for  concrete  targets. 

It  then  follows  that  a rock  penetration  equation  must  be  inherently  less 
accurate  than  a concrete  penetration  equation. 

Due  to  the  limited  penetration  data  base,  the  restrictions  on 
Equation  2 are  as  follows: 

1,  At  best,  the  aceui'acy  of  prediction  is  +20  percent. 

2_.  The  equation  is  valid  only  when  the  calcxdatcd  depth  is 
””  greater  than  three  projectile  diameters. 

3.  For  nearly  intact  x'oek  (RQD  > 90),  the  equation  appears  to  be 
~ applicable  for  projectiles  from  3 to  30  cm  in  diameter. 

t,.  For  rock  with  RQD  < 90,  the  eq\»ation  has  not  been  verified 
” for  projectiles  outside  the  10-  to  30-cm-diam  range. 

If  RQI>  < SO,  the  equation  should  not  be  used  at  all. 

§^,  The  effect  of  nose  shape  seems  to  be  weak,  but  the  equation 
”*  is  not  recommended  for  blunt  or  neai*-blunt  (CRH  <1.5) 
projectiles. 

X*  equation  is  not  valid  if  the  projectile  "muslirooms”  or 
breaks  up. 

8.  The  equation  is  not  valid  if  the  projectile  tumbles,  or  if 
*"  the  penetration  path  is  sharply  eui*vod. 

If  i»aolblo,  the  target  strength  should  be  dotemined  from  static 
unconfined  compi'ession  tests  of  Intact  rock  samples,  and  the  values  used 
for  t and  RCU)  should  oorrospond  to  the  stvmo  borehole. 

Appendix  A contains  information  tliat  cun  bo  used  to  make  rough 
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estimates  of  p , Y and  RQD  in  those  cases  where  only  a word  descrip- 
tion of  the  target  is  available.  However,  in  such  cases  only  very  rough 
estimates  of  penetration  depth  can  he  made.  For  example,  an  uncertainty 
of  +50  percent  in  Y compounded  with  +20  percent  in  RQD  will  produce  an 
error  hand  of  +70  percent  in  the  predicted  depth.  If  there  exist  any 
previous  penetration  data  for  the  site  in  question,  these  data  should  he 
used  as  bench  marks  for  checking  the  results  obtained  from  Equation  2 
and  (indirectly)  the  target  parameters  used  therein. 

Comparison  with  Bomb  Penetration  Data 

Reference  11  (EM  1110-3^*5-^* 3^4)  contains  a rock  penetration  equation 
based  on  data  obtained  by  Livingston  and  Smith  (Reference  12).  In  these 
tests,  inert  l6Q0-pound  AF,  2000-pound  SAP,  and  25,000-pound  SAP  bombs 
were  dropped  on  sandstone  and  weathered  granite  targets  at  impact 
velocities  of  256-361  m/s.  The  bombs  struck  the  targets  at  angles  of 
15-30  degrees  from  the  vertical,  and  the  penetration  paths  were  all 
somewhat  curved.  No  RQD  data  are  available  for  these  tests.  However, 
recorded  descriptions  (and  a few  pictures  of  core  boxes)  suggest  on  RQD 
of  20-lt0  percent  for  the  granite  site  and  50-70  percent  for  the  sand- 
stone site.  Table  3 gives  the  ranges  of  strength  and  density  for  the 
two  target  areas,  and  Table  h presents  the  bomb  chax'acteristlco. 

The  bomb  penetration  depths  were  equated  with  the  total  (curved) 
path  letigthu,  and  the  data  were  nondimensionollsed  (Figure  3)  by  the 
method  used  in  Figure  1.  The  horisontal  error  burs  in  Figure  3 were 
obtained  by  using  the  lover  limits  of  Y and  RQD  to  calculate  the  upper 
bound  of  the  abseissa,  atid  vice  versa.  The  average  density  was  used 
in  each  cose,  and  the  median  values  of  the  abscissa  are  shown  plotted 
with  circles  and  triangles.  The  data  are  badly  scattered  vertically 
(not  due  to  the  nondimensionalisation) , though  Equation  2 is  approxi- 
toately  the  some  as  a least-squar'>a  linear  fit  drawn  tlu*dugh  the  origin. 
The  (vertical)  standard  acviatlon  of  the  data  is  +33  i)ercent.  This 
comparison  demonstrates  both  the  ex{>orimental  uncertainty  and  the  com- 
putational error  that  ai*u  often  met  in  practice. 


Summary 


Equation  2 is  a nondimens ional  empirical  formula  for  calculating 
projectile  penetration  in  massive  rock  deposits.  It  is  expressed  in 
terras  of  standard  engineering  properties,  which  can  he  determined  with- 
out conducting  fvill-scale  penetration  tests.  The  equation  was  derived 
from  a limited  data  base,  and  its  range  of  applicability  should  be  re- 
assessed as  more  penetration  data  become  available.  However,  in  the 
absence  of  a better  prediction  technique.  Equation  2 is  suitable  for 
making  reasonable  estimates  of  final  penetration  depth. 

Recommendat ion 


Equation  2 is  apparently  an  improvement  over  the  rock  penetration 
equation  given  in  Reference  11,  since  (l)  it  fits  the  original  data  on 
which  that  equation  was  based;  (2)  it  fits  more  recent  data;  (3)  it  is 
dimensionally  iiomogeneous ; and  (1()  the  rock  penetrability  is  defined 
in  tex*ffls  of  measui'able  pm'ometers  instead  of  empirical  coefficients. 
Thus,  it  is  recommended  that  the  rock  penetration  equation  (in  any  or 
all  of  its  dimensional,  nondimensional,  or  nomographic  forms)  be  in- 
cluded in  the  nest  revisions  of  References  10  and  11. 


Equation  3 gives  the  dimensional  form  in  metric  u:iits. 
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TABLE  1 PENETRATION  DATA 


Test 

Results 


Target 


Impact 

Velocity 

m/s 


Penetration 

Depth 

cm 


Canfield  and  Clator 

Concrete 

306 

20.3 

tt 

It 

312 

22.9 

tl 

tt 

381 

25.4 

tt 

tt 

453 

36.9 

It 

It 

54i 

41.9 

tt 

It 

602 

59.7 

If 

II 

616 

49.5 

tt 

It 

709 

66.0 

tt 

11 

716 

61.0 

tt 

II 

74i 

69.8 

tt 

II 

773 

73.T 

tt 

II 

809 

, 74.9 

Sandia/DNA 

Welded  Tuff 

372 

222.0 

tt 

It 

4ll 

259.0 

tt 

It 

475 

360.0 

tt 

II 

501 

335.0 

tt 

It 

503 

335.0 

Sandia/DNA 

Sandstone 

444 

357.0 

)i 

•• 

459 

372.0 

Sondia  No.  120-77 

Welded  Agglomerate 

325 

396.0 

Sandia  No.  120-112 

Sandstone 

251 

311.0 

Sondia  No,  120-103 

Sandstone 

268 

305.0' 

Sondia  No.  120-106 

Granite 

262 

381.0 

In  this  tent  the  aaudatone  was  covered  by  76  cm  of  qoH*  so 
the  total  depth  of  penetration  Is  30l  cm.  However,  the  pene- 
tration roslstotico  of  tho  soil  is  negligible  compared  with  that 
of  the  rock. 


11 


■ " " ' 'iffiuatasaSsIfissfaieaB 


O O (M  O t-  CVJ  <M 

O O CO  VO  CO  (O  (O 

H H 


!>» 

+5+5  00 


bO  n 

o 

CVJ 

lA 

OO 

CVJ 

CVJ 

CVJ 

b 0' 

o\ 

O 

0\ 

rH 

r-l 

VO 

a 0) 

« 

« 

t 

• 

* 

♦ 

« 

H Q 

CM 

H 

CVJ 

(H 

CVJ 

CM 

CVJ 

>A 

■fcS. 

,d 

lA 

o 

o 

o 

•6% 

lA 

Cfl 

lA 

. 1 

H 

. 1 

OO 

. 1 

CVJ 

. 1 

H 

. t 

lA 

. • 

H 

u .o  lA  O -:!■  lA  0\  CO  CVJ 


^ O CO  t-  OO  O VO 

CO  VO  CVJ  CVJ  -4- 


to 

y 

<u 

m 

<U 

0) 

«) 

d 

<! 

d 

d 

p 

0 

o 

0 

<u 

'd 

p 

p 

p 

d 

0) 

(U 

u 

•d 

•d 

frt 

tJ 

d 

H 

d 

H 

H 

€ 

0 

a 

a 

o 

ts 

CO 

:* 

CO 

CQ 

o 

O 

o 

lA 

lA 

lA 

lA 

o 

o 

CM 

CM 

CM 

• 

« 

* 

• 

• 

* 

H 

VO 

VO 

<3V 

<3V 

CJv 

CVJ  H H VO  CO  CVJ  W 

VO  lA  lA  CO  00  CO  OO 


(w  VO  - VO 

H H 


CVJ  lA  o O 
CM  CV)  CVJ  CVJ 


p fi 

c^  o o o o o o 


lA  OO  00  O VO  CO  CO 

O O CJV  VO  H H 

CM  CM  CO  ^ VO  VO 


a CO  VO 

q q 

I I T*  1* 

o o o o 

5il  Jil  21  Jii 


«d  «a(  • • • • 

S w o p 5 5 

p Q e ^ to  ia 


d 

1 d 

d 

1 

1 

1 

1 

1 

CJJ 

1 d 

1 CO 

(0 

(0 

•p  p 

d m 

0)  V)  tJ 


P 4J  to 

dd  S 

p 


(0  u (0  0) 
•H  S rl  U 

o 

0)  O P *4 

dl  J4  H P4 

•M  CJ 


« H • 

0)  d ■§ 

MPa 


a a q 

gss  I 


I* 

lf 


TABLE  3 TARGET  PARAMETERS  FOR  BOJffi  PENETRATION 


Material 


Density 

g/cm3 


Strength 

Bars 


Estimated 

RQD 


Sandstone 


2.13  + 0.01 


Weathered 

Granite 


2.73  + 0.02 


3k6  + 112 
630  + 139 


50-70 

20-l;0 


Average  +_  standard  deviation. 


TABLE  k BOMB  PARAMETERS 


Diameter  t cm 


l,600-lb  AP 
2,000-lh  SAP 
25 ♦000-lb  SAP 


748  + 3 
936  + 5 
11,449  i 11 


35.56 

47.65 

58.72 


LEGEND 


SYMBOL 

TARGET 

ROD 

% 

STRENGTH 

BARS 

SOURCE  OF  DATA 

0 

CONCRETE 

100 

345 

CANFIELD  & CLATOR 

L 

WELDED  TUFF 

100 

600 

SANDIA/DNA  TESTS 

□ 

SANDSTONE 

62 

234 

SANDIA/DNA  TESTS 

▲ 

WELDED  AGGLOMERATE 

60 

275 

SANDIA  TEST  NO.  120-77 

• 

SANDSTONE 

37 

489 

SANDIA  TEST  NO.  120-112 

■ 

SANDSTONE 

32 

408 

SANDIA  TEST  NO.  120-103 

? 

GRANITE 

32 

462 

SANDIA  TEST  NO.  120-106 

Figure  1.  I>ata  'baoe  for  Equation  2 
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TABLE  A1  ROCK  QUALITY  DESIGNATION 


RQD,  % 

Rock  Quality 

0-25 

Very  Poor 

25-50 

Poor 

50-75 

Fair 

75-90 

Good 

90-100 

Excellent 

TABLE  A2  ENGINEERING  CLASSIFICATION  FOR  INTACT  ROCK 


Class 

Description 

Compressive 

Strength 

bars 

A 

Very  High  Strength 

Over  2200 

B 

High  Strength 

1100-2200 

C 

Medium  Strength 

550-1100 

D 

Low  Strength 

275-550 

E 

Very  Low  Strength 

Loss  than  275 

18 
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TABLE  A3  COMMON  INTACT  ROCK  DESCRIPTIONS 


Rock  Types 

Typical 

Density 

g/.cm3 

Strength 

Range 

bars 

Soft  Shale 

(clay  shales,  poorly  cemented  silty  or  sandy 
shales) 

2.3 

lU-lliO 

Tuff 

(nonwelded) 

1.9 

14-210 

Sandstone 

(large  grain,  poorly  cemented) 

2.0 

70-210 

Sandstone 

(fine  to  medium  grain) 

2.1 

140-500 

Sandstone 

(very  fine  to  medium  grain,  massive,  well 
cemented) 

2.3 

400-1100 

Shale 

(hard,  tough) 

2.3 

140-800 

Limestone 
(coarse,  porous) 

2.3 

400-850 

Limestone 

(fine  grain,  dense,  massive) 

2,6 

700-14oO 

Basalt 

(vesicular,  glassy) 

2.6 

550-1000 

Basalt 

(massive) 

2.9 

>l4oo 

Quartzite 

2.6 

>l400 

Granite 

(coarse  grain,  altered) 

2.6 

550-1100 

Granite 

(competent,  fine  to  medium  grain) 

2.6 

1000-1900 

Dolomite 

2.5 

700-1400 

■A 
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